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We analyze the sensitivity of /3-decay rates in ^®Ni and iood32gjj ^ correction term in Skyrme 
energy-density functionals (EDF) which modifies the radial shape of the nucleon effective mass. 

This correction is added on top of several Skyrme parametrizations which are selected from their 
effective mass properties and predictions about the stability properties of ^®^Sn. The impact of the 
correction on high-energy collective modes is shown to be moderate. From the comparison of the 
effects induced by the surface-peaked effective mass in the three doubly magic nuclei, it is found 
that is largely impacted by the correction, while ^®Ni and are only moderately affected. 

We conclude that /3-decay rates in these nuclei can be used as a test of different parts of the nuclear 
EDF: ^®Ni and ^°°Sn are mostly sensitive to the particle-hole interaction through the B(GT) values, 
while ^^^^Sn is sensitive to the radial shape of the effective mass. Possible improvements of these 
different parts could therefore be better constrained in the future. 

PACS numbers: 23.40.-s, 26.30.-k, 21.10.Re 


I. INTRODUCTION 

Weak processes such as /3-decay rates, electron cap¬ 
ture, neutrino scattering and absorption play an impor¬ 
tant role during the late evolution of massive stars [l|. 
They are largely responsible for the electron fraction in 
the core during the core-collapse phase and they con¬ 
tinue to plw a determinant role in the nuclear synthesis 
r-process (2|. Because of their great importance in as- 
trophysical applications, weak processes were extensively 
investigated within various approaches. The large-scale 
shell model Monte Carlo (SMMC) method was, for in¬ 
stance, applied to compute /3^ decay rates for stellar 
conditions for more than 100 nuclei in the mass range 
A = 45-65 Si- Recently, mean-field based models have 
been used for the prediction of electron-capture cross 
sections and rates. Finite-temperature charge-exchange 
RPA (CERPA) models based on Skyrme or relativistic 
functionals have been applied to predict electron-capture 
cross sections using different interactions [m. Mean- 
field predictions around the Fermi energy are, however, 
known to suffer from their underestimation of the den¬ 
sity of states. In this work, we explore a small correction 
to the mean field models which increases the density of 
states around the Fermi energy Here, we compare 
the predictions of this model to known experimental val¬ 
ues such as /3 half-lives or collective modes, as a first step 
before using it for astrophysical applications. 

Since the pioneering work of Brown et al. Q it is known 
that the level density around the Fermi energy in sta¬ 
ble nuclei indicates that the in-medium nucleon effective 
mass is close to the bare mass. The description of giant 
resonances such as the giant dipole resonance requires, on 
the other hand, that the nucleon effective mass in the nu¬ 
clear medium should be reduced as compared to its value 
in vacuum 0- Analysis of the momentum dependence 


of the nuclear optical potential also favors an in-medium 
effective mass lower than in vacuum 0[il. 

These apparently diverging properties of the in¬ 
medium effective mass m* can be reconciled by consider¬ 
ing the two contributions to m* : the fc-mass which is also 
called the non-locality mass, and the w-mass which is in¬ 
duced by dynamical correlations such as particle-phonon 
coupling [13l - [l7| . The coupling of the collective modes 
to the single-particle (s.p.) motion is, however, diffi¬ 
cult to perform in a self-consistent approach. One of the 
main problems is coming from the fragmentation of s.p. 
strength which increases exponentially at each iteration 
of the self-consistent method. It has therefore been tried 
to include these correlation directly in the mean field, ei¬ 
ther at the level of the interaction with density-dependent 
gradient terms [l8l| . or, loosing the relation with an in¬ 
teraction, at the level of the nuclear energy density func¬ 
tional (EDF) so as to produce a surface-peaked effective 
mass (SPEM) which, at the same time, does not strongly 
modify the mean-field In this study, we will explore 
the second approach. 

Predictions of /3-decay rates throughout the nuclear 
chart within a consistent microscopic nuclear model are 
difficult. Tuning of models according to the system un¬ 
der study is usually performed, and the description of /3- 
decay rates through a unique microscopic nuclear model 
does not exist. Since /3-decay rates are known to depend 
strongly on the fine structure around the Fermi level, the 
difficulties to have a general description could be related 
to the common issue with mean-field models that the s.p. 
level density around the Fermi level is too low. The in¬ 
crease of the level density, by using for instance a model 
producing a SPEM could, in principle, lead to a bet¬ 
ter description of /3-decay rates throughout the nuclear 
chart. 

In microscopic approaches, calculations of nuclear /3 
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-decay rates are rather complex. Due to phase-space am¬ 
plification effects, the calculated /?-decay rates are sen¬ 
sitive to both nuclear binding energies and /3-strength 
functions. In an appropriate /3-decay model, the correct 
amount of the integral /3-strength should be placed within 
the properly calculated Q^g-window provided that the 
spectral distribution is also close to the ’’true” /3-strength 
function. Furthermore, for consistency the model has 
to yield correct positions and strengths of the Gamow- 
Teller (GT) and first-forbidden resonances in the contin¬ 
uum . Another complication is related with the large- 
scale predictions of nuclear /3-decay rates. Such a pro¬ 
gram is a compromise between accurate results and the 
necessity to cover extended regions of the nuclear chart 
including deformed nuclei or even the region with triple 
prolate-oblate-spherical shape coexistence scenario. In 
this work, we shall consider only the case of spherical nu¬ 
clei. A plausible way to detect a change of the /3-decay 
strength function profile due to higher-order corrections 
could be the analysis of the integral characteristics of /3- 
decay. The half-life is one of such characteristics, being 
sensitive enough to the /3-strength distribution [1^. It 
is worth to analyze first the doubly-magic /3^-unstable 
nuclides, such as since one can use the simpler 

GERPA model. Also, we compare to the most neutron- 
rich {{N—Z)/A = 0.28) doubly-magic nucleus ^®Ni which 
is an important waiting point in the r-process [2^. The 
next step would then be to study the delayed neutron 
and especially delayed multi-neutron emission [^ . This 
is a more difficult task since the delayed neutron emission 
probability (Pn-value) puts an additional constraint 
on the /3-strength distribution in the near-threshold re¬ 
gion. 

This paper is organized as follows. In Sec. II we briefly 
present the modifications to the nuclear EDF which pro¬ 
duce a SPEM, and we describe the protocol used to ad¬ 
just the strength of this correction. In Sec. Ill, we analyze 
the results of the calculations of /3-decay rates in ^®Ni and 
ioo,i32gj^ and the properties of the giant quadrupole res¬ 
onance (GQR) and Gamow-Teller resonance (GTR) of 
^°®Pb. Gonclusions are drawn in Sec. IV. 


II. THE MEAN FIELD MODELS 

Skyrme-type EDF are known to give an accurate de¬ 
scription of masses and charge radii over the whole nu¬ 
clear chart, from Z = 8 up to heavy elements [2^. As 
most of the mean-field approaches, they however lead to 
a s.p. level density around the Fermi surface which is 
lower than the experimental one 0. Here, we introduce 
a correction term to the Skyrme EDF which leads to 
a SPEM and increases the average s.p. level density Q- 
We hereafter present this correction term and then briefly 
describe the calculations of /3-decay rates carried out con¬ 
sistently in the framework of Hartree-Fock-GERPA ap¬ 
proach. 


A. The standard Skyrme functional 

The standard Skyrme functional for the time-even en¬ 
ergy density is expressed as 

'Hsky{r) = X! C?iPo)Pt +C^^Pt^Pt 

+Clptrt + \cUt + CY-^PtV ■ Jt, ( 1 ) 

where the indices t = 0,1 stand for the isoscalar and 
isovector part of the corresponding densities, respec¬ 
tively. For instance, the nucleonic densities po and pi 
are defined as, 

Po{r) = Pn{r)+Pp{r), 

pi{r) = Pn{r)-pp{r), (2) 

where the densities pq {q = n, p) are expressed in terms 
of the s.p. wave functions pf as 

Pq(.r) = (3) 

i 

The kinetic energy and spin-current densities, t* and J*, 
are defined similarly [2^. The coefficients in Eq. (HD 
are constants (see, e.g.. Ref. [2^) except for the coeffi¬ 
cient Cf which depends of the isoscalar density po as: 

Cf (po) = cf (0) + (CHpo.^aO - cf (0)) (, (4) 

\ PO^sat / 

where po^sat is the saturation density in infinite nuclear 
matter. 

The standard Skyrme functional can be separated into 
neutron and proton channels, and neutron and proton 
effective masses are introduced: 

7T7 /TT) 

— = 1 + + Cl)Py + {Cl - CDpq], ( 5 ) 

m* ri^ 

Then, neutron and proton mean fields can be obtained 
(see appendix lAT). 

Among the large number of Skyrme parametrizations, 
we have selected 6 of them based on the following re¬ 
quirements: 

• First, the Skyrme EDF should predict ^®^Sn as a 
,5-unstable nucleus at the mean field level. This is 
based on the common expectation that the Landau 
parameter Gg in the spin-isospin channel is repul¬ 
sive and will shift up the GT strength. 

• Second, we wish to explore different values of ef¬ 
fective mass in the bulk, and different isospin split¬ 
tings of the effective mass. 

In ^®^Sn the first condition can be related to the s.p. 
energy difference between the 7r2(i| and states 

(which contribute mostly to the GT transition towards 



3 


TABLE I: Bulk properties of the selected interactions. 


Skyrme 

P0,sat 

(fm-3) 

Eo 

(MeV) 

Ko 

(MeV) 

Jsym 

(MeV) 

Lsym 

(MeV) 

m* jm 
(MeV) 

Am*/m 

G'o 

SLyIIIO.7 [3^ 

0.153 

-16.33 

361.4 

31.98 

30.78 

0.7 

0.18 

0.30 

SLyIIIO.8 [3^ 

0.153 

-16.32 

368.8 

31.69 

28.24 

0.8 

0.29 

0.33 

SLyIIIO.9 [3^ 

0.153 

-16.31 

374.5 

31.44 

24.75 

0.9 

0.38 

0.34 

/+ [33] 

0.162 

-16.04 

230.0 

32.00 

41.53 

0.7 

0.17 

0.08 

/o [33] 

0.162 

-16.03 

230.0 

32.00 

42.42 

0.7 

0 

-0.01 

/- \m 

0.162 

-16.02 

230.0 

32.00 

43.79 

0.7 

-0.28 

-0.13 


the I’*' state of ^^^Sb). The lowest unperturbed transi¬ 
tion energy is e^ 2 !i| “ ^v 2 d^ ~ where the last 

term stands for the mass difference between the neutron 
and the hydrogen atom, AMn-n = 0.782 MeV. If this 
transition energy is positive at the mean field level - here¬ 
after called the HF transition energy - the system is /3- 
stable since the CERPA correlations could only push it 
up, while it is expected to be actually /3-unstable. Antic¬ 
ipating the discussion of the results in Sec.|TTIlwe observe 
that models having positive HF transition energies pre¬ 
dict /3-decay half-lives which are too large in ^^^Sn. We 
therefore consider only models having a HF energy dif¬ 
ference e^ 2 d^ ~ < 0- This condition is indeed quite 

drastic, and we found that an appreciable number of well 
established Skyrme models do not fulfill it. Among these 
are SHI BSKI4-I7 [H, SKM* [H, SLy4-5 
SKO [ 2 ^. In addition, the models which predict that 
the HF transition energy is larger than 0.782 MeV are: 
RATP [H, SGH [ 23 , LNS [H, LNSI, LNS5 Q, SKIl- 
5 [lOl , SAMi [3l| . These models have therefore not been 
used here. 

For the few remaining models, we restrict our¬ 
selves to the parametrizations SLyHIO.7, SLyHIO.8 and 
SLyHIO.9 im which predict in the bulk nuclear matter 
the effective mass values m*/m=0.7, 0.8 and 0.9, respec¬ 
tively. We have also considered the /_, /q and /+ [s^ 
models which predict an effective mass of 0.7 in sym¬ 
metric matter, with either a positive, zero, or negative 
isospin splitting of the effective mass (ISEM) in neutron 
matter, defined as to*/to — Notice that m*/m 

could be calculated in neutron matter without any am¬ 
biguity: the proton density shall simply be set to zero, 
see for instance Eq. ([5]). 

The bulk properties of the selected interactions are 
given in Table HI It is observed that the saturation den¬ 
sity poj the energy per particle at saturation Eq, and the 
symmetry energy Jsym are very similar for these interac¬ 
tions. The slope of the symmetry energy Lsym is varying 
between 24.75 and 43.79 MeV, which is a rather wide 
range, but these models are still considered as iso-soft 
ones. The incompressibility of SLyHIO.7, SLyHIO.8 and 
SLyHIO.9 is quite large. However, this does not affect 
the processes explored in this work. The main difference 
among these models comes from their effective masses 


and ISEM. The models SLyHIO.7, SLyHIO.8, SLyHIO.9 
have a different effective mass in symmetric matter, and a 
positive ISEM. The models /_, /o and /+ have the same 
effective mass in symmetric matter and different signs for 
the ISEM. The models SLyHIO.7, SLyHIO.8, SLyHIO.9 
and /_|_ have a positive ISEM, as expected from micro¬ 
scopic BHF and DBHF calculations , while /o has no 
splitting and /_ has a negative ISEM. Finally, the val¬ 
ues of the Landau parameter for the models SLyHIO.7, 
SLyHIO.8 and SLyHIO.9 is given by 


G'q = -No 


^^0 + — h) 


( 6 ) 


and for the models /_, /o and /+, 

^to + - t2) 

(7) 

where Nq = 2kpm* is the level density, with kp 
being the Fermi momentum and m* the nucleon effective 
mass. In Eq.(l7]), the parameter is coming from an addi¬ 
tional density-dependent term besides the usual density- 
dependent ta term, and the coefficient in front of the den¬ 
sity dependent terms have been modified w.r.t. the stan¬ 
dard notations . The values of the Landau parameter 
Gq are given in the last column of Table HI At satura¬ 
tion density (p = po), the models SLyHIO.7, SLyHIO.8, 
SLyHIO.9 predict rather large values for Gq w 0.3 — 0.35, 
while the models /_, /q and /+ predict smaller value 
with Gq ~ 0. The forces /_, fp and /+ clearly predict 
not enough positive Gq values [l^]- In addition to the 
different effective masses we therefore expect to observe 
substantial differences between these two sets of models 
in the charge-exchange channel. 


G'o = -No 


B. Surface-peaked effective mass correction 

In Ref. (3, a surface-peaked effective mass correction 
to the Skyrme-type Hamiltonian was proposed with the 
following form. 


An = Go^"^P'>\ {Vpf + {Vpf 


( 8 ) 











4 





r, fm 


r, fm 


r, fm 


FIG. 1: (Color online) Panels (a) and (d): Proton and neutron effective masses in ^°®Pb as functions of radial distance; 
panels (b), (c), (e) and (f): effective-mass difference between the results of the HF calculation with the surface peaked terms 

=-210 MeV fm^° and -420 MeV fm^°), and without. Solid, dashed and dotted lines correspond to HF calculations 
with /o, /- and /+ forces, respectively. 


and the new functional can be written a.s'H = 'Hsky+^'H- 

The first term of Eq.® is designed to modify the effec¬ 
tive mass profile at the nuclear surface, while the second 
term is introduced in order to compensate the effects of 
the first term in the nuclear mean-field. Without the sec¬ 
ond term, the effects of the first term on the mean field 
are too large and drastically limit the possible values for 
the strength of the SPEM, as in Ref. [^. The compensa¬ 
tion was found to be optimal for intermediate mass and 
heavy nuclei if one uses the following constant relation 
between the two new parameters 

CP^iVpf ^ _;^o fm . (9) 

One can expect an impact of the SPEM on the prop¬ 
erties of the lowest quadrupole excitation if the isoscalar 
terms ([5]) are taken into account. On the other hand, the 
energy-weighted sum rule (EWSR) is an integral char¬ 
acteristic and it is particularly sensitive to the giant- 
resonance properties which can be described by the EDE 
without the terms (Ill). In the present work, the values 
of =-210 and -420 MeV fm^° are fixed so that 

the isoscalar quadrupole EWSR in ^°®Pb is modified by 


1% and 2%, respectively. Consequently, a change of less 
than 0.04 of the neutron and proton effective masses at 
the nuclear surface of ^°®Pb is predicted, see Eig. [TJ This 
procedure is slightly different from that used in Ref. 0, 
and it leads to a SPEM less strongly peaked at the sur¬ 
face. We have added the terms dS]) without refitting the 
existing standard parameterizations. Using this pertur¬ 
bative approach, we observe a small change of the binding 
energies which is larger than the tolerance of the proto¬ 
col for the parameter fitting. In particular, in ^°®Pb the 
binding energy changes by 0.35% for the SLyIIIO.9 set, 
0.37% for the SLyIIIO.8 set, 0.38% for the SLyIIIO.7 set 
and 0.45% for the /o, /- and /+ sets. A fine tuning of 
other parameters in order to compensate for these energy 
changes has still to be done. 

In Eig.[T]are shown the effective mass profiles in ^°®Pb 
for the /o, /_ and /+ models where we have considered 
different values of the parameter governing the strength 
of the SPEM, -210 and -420 MeV fm^o. We 

remind that the differences between the models /o, /- 
and /+ are mostly the ISEM in asymmetric matter: /_|_ 
has m* > m* in neutron rich matter, while /_ has 
m* < m*, and /o has m* = m* in the same conditions 
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of isospin asymmetry. The effect of the sign difference of 
the effective mass splitting can also be observed on pan¬ 
els (a) and (d) (without STEM): Since ^°®Pb is a neutron 
rich nucleus, the neutron effective mass is larger than the 
proton one for the f+ model, an opposite effect is found 
for f_, and no effect is observed for /q. Additionally, it 
is observed in panels (b), (c), (e) and (f) that the SPEM 
correction is almost unaffected by the effective mass split¬ 
ting, since the correction is isoscalar. 


C. Calculations of /?-decay rates 

We describe the collective modes in the charge- 
exchange random phase approximation (CERPA) using 
the same Skyrme interactions as above. Making use of 
the finite-rank separable approximation (FRSA) 
for the p-h interaction enables us to perform CERPA 
calculations in very large configuration spaces. Although 
it is well known that the tensor interaction influences 
also the description of the /3“-decay half-lives [s^, in the 
present study the tensor force is neglected in order to 
focus on the impact of the SPEM. 

The experimentally known values of the half-lives put 
an indirect constraint on the calculated GT strength dis¬ 
tributions within the Q^-window. To calculate the half- 
lives an approximation worked out in Ref. (40l | is used. It 
allows one to avoid an implicit calculation of the nuclear 
masses and -values. However, one should realize that 
the related uncertainty in constraining the parent nucleus 
ground state calculated with the chosen Skyrme interac¬ 
tion is transferred to the values of the neutron and proton 
chemical potentials. In the allowed GT approximation, 
the /3^-decay rate is expressed by summing the proba¬ 
bilities of the energetically allowed transitions (in units 
of G^/dTr) weighted with the integrated Fermi function. 
For the /3“-decay case we have: 

y^/3- ^_ £ _ 

{^) ' Efe fo{Z + I, A, E, - E,+ )BiGT )-’ 

( 10 ) 

Ei — E^+ « AMn-H + ~ l-tp ~ (H) 

k 

while for the /3+-decay case this becomes: 


1/2 / \2 ’ 

{^) EkM-Z + hA,E,-E,^)BiGT)t 

( 12 ) 

Ej - E^+ Ki -AMn-H - 2me - Hn + fkp - Ek- (13) 

Here, £)=6147 [4l[ is a constant, GA/Gy=1.25 [4l| is the 
ratio of the weak axial-vector and vector coupling con¬ 
stants, TOe is the positron mass; and are the neu¬ 
tron and proton chemical potentials. Ei is the ground 
state energy of the parent nucleus {Z,A) and E^+ de¬ 
notes a state of the daughter nucleus. The Ek are the I"*" 


TABLE II: Energy differences between the dominant s.p. 
states in ^®^Sn and ^°°Sn. For each Skyrme parameterization, 
the energy difference is calculated with the surface peaked 
term or without = 0). See text for more details. 


Skyrme 

r-rtypf 

'-'0 

(MeV W“) 

“^Sn 

(MeV) 

i°°Sn 

Glgl- 

^Trlgf 

(MeV) 

"«Ni 

G2p|“ 

(MeV) 

SLyIIIO.7 

0 

0.3 

-7.2 

-5.1 

SLyIIIO.7 

-210 

0.2 

-7.3 

-5.2 

SLyIIIO.7 

-420 

0.1 

-7.3 

-5.3 

SLyIIIO.8 

0 

-0.6 

-7.5 

-6.2 

SLyIIIO.8 

-210 

-0.7 

-7.5 

-6.3 

SLyIIIO.8 

-420 

-0.9 

-7.6 

-6.4 

SLyIIIO.9 

0 

-1.3 

-7.7 

-7.0 

SLyIIIO.9 

-210 

-1.4 

-7.7 

-7.2 

SLyIIIO.9 

-420 

-1.6 

-7.8 

-7.3 

/+ 

0 

-0.6 

-5.9 

-5.8 

/+ 

-210 

-0.7 

-6.0 

-5.9 

/+ 

-420 

-0.8 

-6.2 

-6.0 

/o 

0 

-0.5 

-6.0 

-5.6 

fo 

-210 

-0.6 

-6.1 

-5.7 

/o 

-420 

-0.7 

-6.2 

-5.8 

/- 

0 

-0.4 

-6.2 

-5.4 

/- 

-210 

-0.5 

-6.3 

-5.5 

/- 

-420 

-0.6 

-6.5 

-5.6 


eigenvalues of the GERPA equations. The GERPA wave 
functions allow us to determine GT transitions whose op¬ 
erator is d±=Y,im 


B{GT)t= {N±l,ZTlAt\0£N,Z-Q%) 


(14) 


Expressions (Unil-dlll) will be used in the next section 
to calculate the /3-decay rates and the collective modes. 
All the calculations are performed without any quenching 
factor. 


III. RESULTS FOR COLLECTIVE MODES AND 
/3-DECAY RATES 

We now analyze first the results of the /3-decay rates 
which are sensitive to the low-energy part of the CERPA 
strength, and then the GT collective modes. The effects 
of the SPEM will be discussed. 

The p-h interaction in the spin-isospin channel is as¬ 
sumed of the following form: 

U(ri,r2) = Alo'^Go(r'i)(T(^)-(T(^V(^)-T^^^(5(ri-r2) (15) 

where and are the spin and isospin operators. 
As expected, the largest contribution to the calculated 
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TABLE III: SPEM effects on P -decay properties of 
Data are from Ref. [4^. 


Skyrme 

n-riypp 

(MeV fmi“) 

Ei — 

(MeV) 

i3(GT)r 

T\/2 

(s) 

SLyIIIO.7 

0 

0.07 

2.6 

389400 

SLyIIIO.7 

-210 

0.21 

2.7 

9840 

SLyIIIO.7 

-420 

0.34 

2.7 

1930 

SLyIIIO.8 

0 

0.97 

2.5 

57 

SLyIIIO.8 

-210 

1.11 

2.5 

33 

SLyIIIO.8 

-420 

1.26 

2.6 

21 

SLyIIIO.9 

0 

1.70 

2.4 

6.7 

SLyIIIO.9 

-210 

1.84 

2.5 

4.7 

SLyIIIO.9 

-420 

2.01 

2.6 

3.3 

/+ 

0 

1.12 

4.6 

18 

/+ 

-210 

1.25 

4.6 

12 

/+ 

-420 

1.36 

4.6 

8.5 

fo 

0 

1.14 

5.9 

13 

/o 

-210 

1.27 

5.8 

00 

bo 

fo 

-420 

1.37 

5.8 

6.5 

f- 

0 

1.12 

8.8 

6.4 

f- 

-210 

1.25 

8.7 

5.0 

f- 

-420 

1.36 

8.6 

3.6 

Expt. 


1.794T0.009 


39.7T0.8 


TABLE IV: SPEM effects on P -decay properties of ’^^Ni. 
Data are from Refs. [43|. 


Skyrme 

(MeV fm“) 

Ei — E,+ 
-^1 

(MeV) 

B{GT)p 

Ti/2 

(s) 

SLyIIIO.7 

0 

5.49 

1.0 

0.157 

SLyIIIO.7 

-210 

5.61 

1.1 

0.140 

SLyIIIO.7 

-420 

5.74 

1.1 

0.121 

SLyIIIO.8 

0 

6.60 

1.0 

0.057 

SLyIIIO.8 

-210 

6.73 

1.0 

0.051 

SLyIIIO.8 

-420 

6.87 

1.0 

0.045 

SLyIIIO.9 

0 

7.48 

1.0 

0.025 

SLyIIIO.9 

-210 

7.61 

1.0 

0.023 

SLyIIIO.9 

-420 

7.79 

1.0 

0.020 

/+ 

0 

6.40 

1.9 

0.031 

/+ 

-210 

6.51 

1.9 

0.028 

/+ 

-420 

6.60 

1.9 

0.027 

fo 

0 

6.33 

2.6 

0.020 

fo 

-210 

6.44 

2.5 

0.019 

fo 

-420 

6.52 

2.5 

0.018 

f- 

0 

6.33 

3.9 

0.010 

f- 

-210 

6.42 

3.9 

0.009 

f- 

-420 

6.50 

3.8 

0.009 

Expt. 



0.1222T0.0051 


/3^-decay half-life comes from the I]'' state, the struc¬ 
ture of which is dominated by one unperturbed configu¬ 
ration. They are the Ip-lh configurations {7r2d|, 

7 rl(7|} and {7r2p|, v2p^} of ^^^Sn, ^™Sn and "^^Ni, 
respectively. In other words, the 1^ state is non-collective 
and, therefore, the /3-decay is related to the lowest unper¬ 
turbed 1^ energy. We first examine the s.p. energy dif¬ 
ferences given in Table HIl for the selected Skyrme models 
and for various strength of the SPEM parameter Cg. 
They are small (about 1 MeV) in ^^^Sn but rather large in 
^®Ni and ^°°Sn (5 to 8 MeV). In ^°°Sn, the energy differ¬ 
ences without SPEM are mostly sensitive to the Coulomb 
component of the EDF, with a small additional effect 
due to the effective mass (the larger effective mass, the 
smaller the energy difference). In ^^^Sn, the energy differ¬ 
ence is related mostly to the symmetry energy: the larger 
the symmetry energy, going from SLyIIIO.9 to SLyIIIO.7 
for instance, the larger the energy difference. In addition, 
the increase of the effective mass also contributes, with 
a smaller impact, to the decrease the energy difference, 
as it can be deduced from the comparison of the energy 
difference for the forces /_, /o, and which has increas¬ 
ing effective mass in neutron rich matter, see Fig. [TJ It 
can be seen (cf. Table |lT| that the shifts in the energy 
differences between the cases without, and with maximal 
SPEM (Cq*-^^^ =-420 MeV fm^°) are almost constant 
and independent of the models considered. It varies by 
about 0.3 MeV in ^^^Sn and in ^®Ni. From Table II, we 
can anticipate that the SPEM will have a larger impact 
on the calculation of the /3 half-life of ^^^Sn and a weaker 
one in the case of ^®Ni and ^™Sn. In ^^^Sn, the experi¬ 
mental value is -1.305 MeV [d^. No data exist for ^®Ni 
and i°°Sn. 

The Ei — E^+ energies, the B{GT)^ values and P~- 
decay half-lives of ^^^Sn and ^®Ni are given in Tables Hill 
and lIVi the ,5“''-decay properties of ^^'^Sn in Table IVl The 
evolution of the transition energies and the B{GT)^ val¬ 
ues is reflected in the half-life behaviour, see Eqs. m 
and ((HI). As in Table [TTl the results shown in Ta¬ 
bles mil HV] and |V] correspond to the selected interac¬ 
tions with and without the SPEM represented by the 

value of the parameter . In the case of ^^^Sn, 

the model SLyIIIO.7 predicts positive energy differences 
for the dominant transition of the /3-decay half-lives (see 
Table HIl) ! and it leads to half-lives which are much larger 
than the experimental value, as anticipated. 

One can see from Tables Hill llVI andIVlthat the /3-decay 
half-lives are much more sensitive to the effective mass 
distribution in the case of the low-Q^g nucleus ^^^Sn than 
in ^™Sn and ^®Ni. For ^^^Sn, a strong decrease of the 
half-life can be directly correlated to either the increase of 
the effective mass in symmetric matter, or to the increase 
of the SPEM, while in ^®Ni and ^°°Sn, the correlation, 
while still being present, is much less pronounced. This 
can be easily understood from the energy difference of the 
most important transition given in Table HIl The energy 
differences are much smaller in the case of ^®^Sn than 












TABLE V: SPEM effects on /^"'"-decay properties of 
Data are from Ref. [4^. 


Skyrme 

^t(Vp)2 

(MeV W°) 

Ei — E.+ 

(MeV) 

B{GT)+ 

Ti/2 

(s) 

SLyIIIO.7 

0 

4.26 

15.2 

0.232 

SLyIIIO.7 

-210 

4.38 

15.2 

0.221 

SLyIIIO.7 

-420 

4.42 

15.2 

0.213 

SLyIIIO.8 

0 

4.60 

15.1 

0.178 

SLyIIIO.8 

-210 

4.64 

15.1 

0.172 

SLyIIIO.8 

-420 

4.67 

15.0 

0.167 

SLyIIIO.9 

0 

4.86 

15.1 

0.138 

SLyIIIO.9 

-210 

4.90 

15.0 

0.134 

SLyIIIO.9 

-420 

4.92 

15.0 

0.131 

/+ 

0 

3.47 

16.3 

0.593 

/+ 

-210 

3.62 

16.3 

0.492 

/+ 

-420 

3.72 

16.2 

0.433 

fo 

0 

3.80 

16.8 

0.381 

fo 

-210 

3.94 

16.8 

0.323 

fo 

-420 

4.04 

16.7 

0.290 

f- 

0 

4.38 

17.6 

0.190 

f- 

-210 

4.51 

17.5 

0.168 

f- 

-420 

4.60 

17.5 

0.154 

Expt. 


3.08T0.34 


1.16T0.20 


TABLE VI: SPEM effects on the energy and B(E2)-value for 
the up-transition to the first 2^ state in ^'^®Pb. Data are from 
Ref. ra. 


Skyrme 

^0 

(MeV fmi°) 

Energy 

(MeV) 

B(A2;0+ ^2+) 
(e^fm'^) 

/+ 

0 

5.12 

3130 

/+ 

-210 

5.09 

2530 

/+ 

-420 

5.09 

2180 

fo 

0 

5.13 

3250 

fo 

-210 

5.09 

2650 

fo 

-420 

5.09 

2310 

f- 

0 

5.09 

3440 

f- 

-210 

5.06 

2850 

f- 

-420 

5.06 

2500 

Expt. 


4.09 

3000T300 


in the case of and ^®Ni, which makes the /3-decay 
half-lives more sensitive to a small modihcation of the 
s.p. energies induced by the SPEM. 

Let us examine whether the SPEM could improve the 
agreement between the model predictions and the exper¬ 
imental values. As one can see from Tables III, IV and 
V the inclusion of the terms (0) leads to minor effects 
on the B{GT)^ values. We hnd that the SPEM induces 
an increase of the transition energies and it results in a 
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FIG. 2: (Color online) The quadrupole strength distribution 
of ^°®Pb. Solid, dashed and dotted lines correspond to RPA 
calculations with /o, /_ and /+ models, respectively. The 
experimental centroid of the GQR is at 10.89T0.30 MeV 


decrease of the half-lives. We first concentrate on the 
models SLyIIIO.7, SLyIIIO.8 and SLyIIIO.9, which corre¬ 
spond to different values of the effective mass in sym¬ 
metric matter. From the comparison of the theoretical 
predictions with the experimental half-lives shown in Ta¬ 
bles |nTl |TV] and El it is difficult to conclude which model 
is better: For ^^^Sn, the model SLyIIIO.8 is preferred, for 
78Ni and i™Sn, it is SLyIIIO.7. Now, if we concentrate on 
the models f+, fo and f_, it is f+ which always comes the 
closest to the experimental value. This indicates that, in 
addition to the effective mass, the residual interaction is 
very important. It was already anticipated that the value 
of the Landau parameter Gg for the selected models (cf. 
Table m could have an impact on charge-exchange related 
observables. For the models f+, fo and f_ the values of 
Gq are too small. Since the impact of the SPEM on the 
/3-decay rates in ^®Ni and ^°°Sn is quite small, these nu¬ 
clei could be, in the future, used to calibrate the residual 
interaction almost independently from the profile of the 
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FIG. 3: (Color online) Same as Fig. [21 for the GT strength 
distribution obtained within the CERPA. The experimental 
centroids of the GTR is at 19.2 MeV [i?!. 1^. 



effective mass. The modification of Gg could be obtained 
from a refitting of the Skyrme functional with a different 
value of the strength of the SPEM. One could increase 
Gg by about 0.1-0.2 by introducing the spin-density de¬ 
pendent extension of the Skyrme model (d^ . IHoj ]. This 
will be left for future investigations. 

Up to this point, we have mostly focused on the rela¬ 
tion between the SPEM and the low energy part of the 
strength, since it represents the main contribution to the 
/3-decay rates. We now turn to the higher energy part 
and show in Figs. [2] and [3] the effects of the SPEM on 
the properties of the GQR and GTR in ^°®Pb. In the 
figures, the calculated strength distributions are folded 
out with a Lorentzian distribution of 1 MeV width. The 
excitation energies refer to the ground state of the parent 
nucleus ^°®Pb. The arrows indicate the maxima of the 
strength distributions corresponding to the case of the 
/+ model and Gg^^^^ =0 MeV fm^°. Since the isoscalar 
quadrupole EWSR is changed by only about 1% by the 
SPEM, we expect the collective modes at higher energy 


to be only marginally impacted. 

The GQR strength distribution consists mostly of a 
main peak. Comparing the cases without SPEM and 
with maximal SPEM, we find that the peak is shifted 
down by about 500 keV. One can notice that the GQR 
strength distribution is almost identical for the three 
models /g, /+ and /_. As can be seen from Ta¬ 
ble IVIl while the 2^ energy is practically unaffected by 
the SPEM, the B(E2) value is decreasing as the SPEM 
increases. Some overestimate of the experimental en¬ 
ergy indicates that there is room for the two-phonon ef¬ 
fects [5l[. 

The GTR is much more fragmented than the GQR, 
as seen in Fig. [3] The strength distribution is globally 
shifted up as the isospin splitting is going from positive 
(f+) to negative values (f_). As in the case of the GQR, 
the high energy peaks of the strength distribution are 
shifted to lower energies (by about 500 keV) as the SPEM 
gets larger. This is an effect of the slight increase of the 
level density induced by the SPEM. 


IV. CONCLUSIONS 

Starting from different Skyrme EDF which predict 
^^^Sn /3-unstable, we have studied the effects of intro¬ 
ducing a surface-peaked effective mass on top of existing 
Skyrme models. The main effect of this additional term 
is a compression of the s.p. level spacing around the 
Fermi level, or equivalently, an increase of the level den¬ 
sity. This systematically increases the /3-decay rates (i.e., 
decreases the half-lives). The collective modes at higher 
energy are only slightly impacted by the SPEM. 

This work is a first step towards improving Skyrme 
functionals by adding extra terms to the energy func¬ 
tional. Our motivation is based on both having a better 
agreement with nuclear data, and also predicting weak 
transition rates for astrophysical applications. The re¬ 
sults of our analysis allow for a better understanding of 
the effects at play. The /3-decay rates in doubly-magic 
unstable nuclei ™Ni) are indeed very sensi¬ 

tive both to the s.p. energies and residual interactions, 
and none of the Skyrme models selected in this work are 
fully satisfactory in this respect. From our analysis, we 
have however identified two nuclei (^™Sn, ^®Ni) where 
the /3-decay half-lives are only weakly impacted by the 
SPEM. They can be considered as good benchmark nu¬ 
clei since they potentially offer the possibility to calibrate 
the residual interaction, with a weak influence of the ef¬ 
fective mass. In a complementary approach, ^®^Sn could 
be used to test different strengths of the SPEM, for a 
fixed residual interaction. 

The tensor force has not been considered in this work, 
although it can affect the neutron-proton s.p. energies 
in some cases. We have aimed at understanding just the 
contribution of the SPEM to the ,0-decay and GT mode 
in order to disentangle the respective roles of the effec¬ 
tive mass and the residual interaction. An additional 
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modification of the Skyrme functional was proposed ear¬ 
lier in order to stabilize the nuclear matter equation of 
state [i^, [13 ■ It has been recently used in nuclei and, 
since it brings an additional repulsive term to the Gq 
Landau parameter, it was shown to shift the centroids 
of the GT collective mode to higher energies by a few 
hundreds keV up to 1 MeV [s^. In the future, we plan 
to explore the predictions of a general mean held model 
including all these ingredients, and to compare them to 
known experimental data, as done in this work. These 
calibration processes are important to set-up boundaries 
for the additional parameters before making predictions 
for astrophysical cases. 
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Appendix A: Decomposition of the Skyrme 
functional into neutron and proton channels 


where the central potential is given by 

Pf y (r) = VP{r) + p 7 (r) + P/(r). (A7) 

Here, the central-density potential is given by: 

= (Go" + G[)r, + (CJ-Cr)r,- 
+ 2[(C'o -I- Ci)pq + (Cg - Ci)pq] 

+ ^^(C'o + C'l )P9 + “ C'fjpqPg], 

(A8) 

the central-gradient potential by: 

P7(r) = 2(Co^" + + 2(Co^" - 

(A9) 

and the central-J potential by: 

Vq^{r) = • Jq + (Cg^^ - • Jq. 

(AlO) 

The spin-orbit potential is: 


Here, the Skyrme functional is expressed in terms of 
the neutron and proton densities instead of the isoscalar 
and isovector densities, 

nskv{.r)= + (Al) 

q—n,p 

where the different terms of the energy density are 

K = ^fq'^^rq + {C^,+C^,)pl + {C^,-Cl)pqPq 

(A2) 

h^q = -(Cg^'’ + GfO(Vp,f-(C'g^'’-Cf'’)Vp,•Vp,- 

(A3) 

H = \{C^o+C()Jl + \[C^,-C()JqJq 

- [(Cg^-' + CY^)Vpq + - CY^)yPq] ' Jq, 

(A4) 

and the effective mass factor = m/m* is defined as 

/f" = 1 + TP [(^0 + Cl)Pq + (Q - Cl)pq] . (A5) 

By functional derivation the one-body Hamiltonian T-Lq 
is obtained as, 

-Hq = -|^V-/fy(r)V + P,(r)- 

cr' 

(A6) 


Wq{r) = + CY^)ypq-{C^^ -CY'’)ypq 

+ (Co^ + C()Jq + (Cg^ - C()Jq. (All) 


Appendix B: Modification of the mean-field 
equations induced by the SPEM 

The kinetic energy correction induced by the effective 
mass in Eq. (IA2I) is now given by fq = where 

/corr ^ (Vp(r))^ , (Bl) 

and the mean field central potential (TATI) reads: 

= + , (B2) 

where I/j®*'y(r) is the mean field deduced from the Skyrme 
interaction, e.g. Eq. (TATI) . and P'^°’'’'(r) is the correction 
term induced by Eq. (jS)): 

^corr(r) ^ (T(r)vV(r) -f VT(r)Vp(r)) 

-2C^ (^p(r)(Vp(r))2-f p(r)2vV(r)) . 

(B3) 
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